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Abstract 

A transformation is devised to convert any lattice Dirac fermion op- 
erator into a Ginsparg- Wilson Dirac fermion operator. For the standard 
Wilson-Dirac lattice fermion operator, the transformed new operator is 
local, free of 0{a) lattice artifacts, has correct axial anomaly in the 
trivial sector, and is not plagued by the notorious problems ( e.g., addi- 
tive mass renormalization ) which occur to the standard Wilson-Dirac 
lattice fermion operator. 
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1 Introduction 

To formulate chiral fermions on the lattice, one must take into account the 
constraints imposed by the Nielson-Ninomiya no-go theorem J5J. It asserts 
that any lattice Dirac operator D in the free fermion limit must violate at 
least one of the four basic properties [(i)-(iv) listed below ] of massless Dirac 
fermion in continuum. Nevertheless, we can incorporate the gauge interactions 
in the no-go theorem as follows. Since the gauge link variables are trivial in 
the Dirac space, the chiral symmetry of D is not affected by turning on a 
background gauge field. In other words, if {Z),7s} is zero in the free fermion 
limit, then it remains to be zero even when the gauge links are different from 
the identity. Moreover, the presence of gauge interactions does not improve the 
locality of D. Therefore, we can assert that any gauge covariant Dirac operator 
Dona finite lattice must violate at least one of the following properties : 

(i) D is chirally symmetric. 



(ii) D is local. 

( \D(x, y)\ ~ exp(— \x — y\/l) with I ~ a; or D(x, y) = for \x — y\ > z, 
where z is much less than the size of the lattice. ) 

(iii) In the free fermion limit, D is free of species doublings. 

( The free fermion propagator D~ 1 (p) has only one simple pole at the 
origin p = in the Brillouin zone. ) 

(iv) In the free fermion limit, D has correct continuum behavior. 
( In the limit a — > 0, D(p) ~ ijnP^ around p — 0. ) 

During the last two years, it has become clear that the proper way to 
circumvent the no-go theorem is to break the chiral symmetry of D (i) at 
finite lattice spacing, according to the Ginsparg- Wilson relation 



where R is any positive definite Hermitian operator which is local in the po- 
sition space and commutes with 75. However, the GW relation should be 
regarded as a generalized chiral symmetry which contains the usual chiral 
symmetry (i) in the continuum limit a — > 0. 

Since the massless Dirac operator in continuum is chirally symmetric ( 
1^75 + 75P = ) and antihermitian ( T>^ = —T> ), so it is 75-Hermitian 
( X>t = 75D75 ). Thus, we require that the lattice Dirac operator D also 
preserves this symmetry at any lattice spacing, i.e., 



(£75 + 75,0 = 0. ) 



1)75(1 - aRD) + (1 - aDR)j 5 D = , 



(1) 



75^75 • 



(2) 
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The first explicit and physical D satisfying the GW relation ([I]) ( with R = 
1/2 ) and 75-Hermiticity is the Overlap Dirac operator D = a _1 (I + 7se(-ff w )) 
which was obtained by Neuberger [0] in the framework of the Overlap formal- 
ism [|, |[, before the GW relation was rediscovered. Since the Overlap was 
essentially motivated by the Domain- Wall fermion |J which was originally 
apart from the Ginsparg- Wilson fermion, this seems to indicate that the Over- 
lap plays a more fundamental role in resolving the problem of chiral fermions 
on the lattice than the GW relation. Thus the Overlap possesses the general- 
ized chiral symmetry (|l|) on the lattice as one of its ingredients. This in turn 
suggests that any nonperturbative formulation of chiral fermions must comply 
with the exact chiral symmetry, the GW relation ([!]). 

One of the salient features of any GW Dirac operator satisfying (ii)-(iv) 
and (H) is that it is not plagued by the notorious problems [|7| ( e.g., additive 
mass renormalization, etc. ) which usually occur to the Wilson-Dirac lattice 
fermion operator. However, if one uses the Overlap Dirac operator in lattice 
QCD, one may encounter the technical problem of taking the square root of a 
huge positive definite Hermitian matrix H^, which would become quite time- 
consuming when the condition number of gets too large. Therefore, both 
in principle and in practice, it is important to understand to what extent one 
can construct a topologically proper D which satisfies all physical constraints 
(ii)-(iv), ® and (|). 

The general solution of the Ginsparg- Wilson relation (|l|) can be formally 
written as J9| |TU| 



D = D c (l + aRD c )~ = (1 + aD c R)~ D c (3) 

where D c is any chirally symmetric ( -D c 7s + 75-Dc = ) Dirac operator which 
must violate at least one of the three properties (ii)-(iv) listed above. Now 
we require D c to satisfy (iii) and (iv), but violate (ii) ( i.e., D c is nonlocal ), 
since (W) can transform the nonlocal D r into a local Dona finite lattice for 



R = rl with r in the proper range [11], and also preserves the properties (iii) 
and (iv). 

Then the 75-Hermiticity of D (|2|) is equivalent to the 75-Hermiticity of D c 

Dl = 75^c75 , (4) 

due to the relation 

D c = D(l - aRD)- 1 = (1 - aDR^D , (5) 

which is the inverse transformation of (|3[). Then the chiral symmetry and the 
75-Hermiticity of D c together implies that D c is antihermitian 

D\ = ~D C , (6) 
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which is in agreement with the antihermiticity of the massless Dirac operator 
in continuum. Then there exists one to one correspondence between D c and a 
unitary operator V such that 

D c = Ma- 1 (l + V)(l-V)-\ (7) 
V = (aD c - M)(aD c + M)-\ (8) 

where M is a parameter. Then V also satisfies the 75-Hermiticity, = 75^75. 
Substituting (^) into (|), we obtain || 

D = MaT l {l + V)[(l -V) + MR(1 + V)}' 1 . (9) 

Note that D is well-defined even if D c has singularities in nontrivial gauge 
backgrounds. In this case, the index of D c is defined to be the index of (I+V), 
the numerator on the r.h.s. of (0). Evidently, the zero modes and the index of 
D are invariant for any R in (|3]) or (|9]) || [H|. That is, a zero mode of D c is 
also a zero mode of D and vice versa, hence, 

n + {D c ) =n+(D), n^{D c ) = n_(D), (10) 

index(D c ) = n_(D c ) - n + (D c ) = n^{D) - n+{D) = index(D) , (11) 

where n + (nJ) denotes the number of zero modes of +1(— 1) chirality. 

The conditions (ii)-(iv), (|l|) and (0) constitute the neccesary conditions 
for D to reproduce the continuum axial anomaly on a finite lattice. For the 
trivial sector, they are sufficient to guarantee that the correct axial anomaly 
can be recovered on the lattice. 

However, for topologically nontrivial gauge backgrounds, the index as well 
as the axial anomaly of a lattice Dirac operator D depends on its topolog- 
ical characteristics If D is topologically proper ( i.e., satisfying the 
Atiyah-Singer index theorem for any gauge background satisfying the topo- 
logical bound [TI| ), then the sum of the axial anomaly of D over all sites on a 
finite lattice is equal to the topological charge of the gauge background. Then 
it follows that the axial anomaly of D would agree with the topological charge 
density of the gauge background if D is local. 

Now the central problem is how to construct a chirally symmetric and non- 
local D c which satisfies (i), (iii), (iv), and (01). Furthermore we also require 
that D c is topologically proper. These constitute the necessary requirements 
(J for D c to enter ([3|) such that D could provide a nonperturbative regu- 



larization for a massless Dirac fermion interacting with a background gauge 
field. A systematic construction of D c satisfying (i), (iii), (iv), and ([|) has 
been discussed in ref. [[13]. However, a general prescription for constructing a 
topologically proper D c ( D ) still remains an unsolved problem. 

Nevertheless, given any lattice Dirac operator D, one can always transform 
D into a GW Dirac operator D' which preserves all essential properties of D. 
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In other words, if D is topologically proper and satisfies (ii)-(iv) and (0), then 
the GW Dirac operator D' also possesses these properties. 

The outline of this paper is as follows. In Section 2, we define the topolog- 
ically invariant transformations for lattice Dirac operators. The transfor- 
mation which can convert any lattice Dirac operator into a GW Dirac operator 
is derived. In Section 3, we apply this transformation to the Wilson-Dirac lat- 
tice fermion operator and obtain a new GW Dirac operator which is local, free 
of 0(a) lattice artifacts and has correct axial anomaly in the trivial sector. In 
Section 4, we discuss and summarize. 



2 Topologically invariant transformation 

Given any lattice Dirac operator D, in general, there are many different ways 
to extract its chirally symmetric part. For example, one can construct 

D s = ^(D- l5 D l5 ) (12) 



which is chirally symmetric ( -D s 7 5 + 7s-D s = ). However, ([12]) does not 
necessarily preserve the property (iii). For example, if one applies this trans- 



formation to the Wilson-Dirac operator Dw = 7^ + W (20), one obtains 
D s = 7^/x, the naive fermion operator which suffers from the species dou- 
blings. Although Dw is free of species doublings in the continuum limit, the 
transformation ( Tl"2|) cannot preserve this property since D s satisfies (i), (ii) 
and (iv), thus it must violate (iii) as a consequence of the no-go theorem. 
Therefore, we need a transformaton which preserves the properties (iii), (iv), 
(0) and (0), but exchanges the locality of D for its chiral symmetry at finite 
lattice spacing. 

Consider the topologically invariant transformation |T5[] on any 75-Hermitian 

D, 

T(R): D^D ! = T(R)[D] = D(l + aRD)- l = (l + aDR)- 1 D (13) 

where R is a Hermitian operator which commutes with 75 and satisfies 
det(l + aRD) 7^ 0. Then D' is also 75-Hermitian. 

It is obvious that a zero mode of D is a zero mode of D' and vice versa. 
Suppose is a zero mode of D, i.e., Dip = 0, then D'<p = (1 + aDR)~ l D<p = 0. 
On the other hand, if is a zero mode of D', i.e., D'cf) = 0, then we can use 
the relation (1 + aDR)D' = D to obtain Dip = 0. Therefore 

= <J=> D'S = . 



Moreover, if the zero modes of D have definite chirality, i.e., Dip = implies 
-^750 = 0, then the chirality of <p is preserved under the transformation (|13|), 
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i.e., -D'750 = 0. In this case, we have 

D(f) ± = D'(j)± = , 

n ± {D>) = n ± (D) , 
index (D 1 ) = index (D) . 

In the trivial sector, D is nonsingular except possibly some "exceptional" 
configurations of zero measure, thus ([13]) is equivalent to 



T{R) : D -»• D' such that D'^ 1 = D 1 + aR . (14) 

Formally, there exists a set of transformations {T(R)} which form an 
abelian group with group multiplication 

T(R 1 ) o T(R 2 ) [D] = T{R±) [ T(R 2 )[D] ] = T(R 1 + R 2 )[D] , (15) 



where the last equality follows immediately from (0). Explicitly, one checks : 

(a) the closure property is satisfied since the sum of any two (75-)Hermitian 
operators is (7 5 -)Hermitian; 

(b) the identity element is T(0); 

(c) the inverse of T(R) is T(—R); 

(d) the associative law 

(T(i? 1 ) o T(R 2 )) o T(R 3 ) = T(R 1 ) o (T(R 2 ) o T(R 3 )) 
= T(R 1 + R 2 + R 3 ) 

is also satisfied. 

The chiral limit of the transformation (0) is at 



R = + ^D- 1 ^) = R c (16) 

which gives the chirally symmetric D c 

D c = T(R C )[D] = 2 l5 D( l5 D - D^D . (17) 

Substituting D c (|17|) into (|]), we obtain a GW Dirac operator 

D' = 2 l5 D( l5 D - D l5 + 2aDR l5 D)- 1 D , (18) 

which satisfies the GW relation 

D' l5 + 75 /r = 2aD'R l5 D' , 

where R is a positive definite Hermitian operator which is local in the position 
space, commutes with 75, and is chosen such that the inverse operator on the 
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r.h.s. of fll8"|) is well-defined. In the following, we shall denote the GW Dirac 
operator ([18]) by 



D' = T{R) [D B ] = T{R) o T{R C ) [D] = T{R + R c ) [D] . 

Note that if D has zero modes in topologically nontrivial gauge back- 
grounds, then D c ( |17|) has singularities and is not well-defined, however, the 
GW Dirac operator D' ([Tj]) is still well-defined. For example, consider the 
75-hermitian Dirac operator D = a _1 (l + V), where V is a unitary operator 
having real eigenvalues ±1 in nontrivial gauge backgrounds. Then it can be 
shown that the real ( ±1 ) eigenmodes of V have definite chirality and sat- 
isfy the chirality sum rule, while each complex eigenmode has zero chirality. 
Thus each +1 eigenmode must be accompanied by a —1 eigenmode of opposite 
chirality, and vice versa. Now Eq. (|17|) gives 



which has a pole if D has a zero mode. However, Eq. (|T^) gives a well-defined 
GW Dirac operator 

D' = 2a- 1 (l + V)[(l -V) + 2R(1 + V)}' 1 . 

It is evident that for any two lattice Dirac operators and satisfying 
(0), their corresponding chiral limits obtained from (0), say, D)p and Df\ 
are in general different. However, they are related by the transformation 

DV = Y.TiD?T} (19) 

i 

where each Tj commutes with 75. In general, given any two T>« and Df \ it 
is nontrivial to obtain all T in fll9l). 



3 Chirally invariant Wilson-Dirac operator 

The Wilson-Dirac lattice fermion operator || can be written as 

D w = 7/A + W (20) 

where 

t ^ x ' ^ = ~^ u ^ 5 *+^y ~ U t(y) s x-»J > ( 21 ) 
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a^al + a v o\ = 25^ , (23) 



and W is the Wilson term 

W(x, y) = ^- E [ 2( ^ - W<W,, - U lW*-i>,y] ■ ( 24 ) 

The color and Dirac indices have been suppressed in (|20|). The first term on 
the r.h.s. of (|20D is the naive fermion operator which satisfies properties (i), 
(ii) and (iv) but violates (iii) since it has 2 d — 1 fermion doubled modes. The 
purpose of the Wilson term is to give each doubled mode a mass of ~ 1/a 
such that in the continuum limit each doubled mode becomes infinitely heavy 
and decouples from the fermion propagator. However, the introduction of 
the Wilson term has serious drawbacks. It causes 0(a) artifacts and also 
leads to the notorious problems such as additive fermion mass renormalization, 
vector current renormalization, and mixings between operators in different 
chiral representations. 

The free fermion propagator of Dw in momentum space is 

t 2 w 

where = ia~ l sin(p M a), t 2 = a~ 2 J2u sin 2 (p M a) and w = 2a~ l Y^ fJ ,sm 2 (p fi a/2). 
We note that in the first term of (^5)), the doubled modes are decoupled due to 
the vanishing of the factor t 2 / iw 2 +t 2 ) at the 2 d — 1 corners of the Brillouin zone. 
However, the second term in (25) breaks the chiral symmetry explicitly, and 



is the source of additive mass renormalization and other notorious problems. 
Any satisfactory solution to all these problems must get rid of the second term 
entirely, while keeping the first term intact. 

Applying the topologically invariant transformation ([17]) to Dw, we obtain 
the chirally symmetric D c , 

D c = - W (7^) _1 W (26) 

which is antihermitian, nonlocal ( due to the second term ), free of 0(a) lattice 
artifacts, and satisfies (iii), (iv) and (f|). 

The free fermion propagator of ( p6|) in momentum space is 

D-\p) = (^t,)- 1 ^— (27) 

which is exactly the first term of (f25|). So, we have "chiraled away" the second 
term of (|25| ) through the transformation (p!7|). 

Substituting fl26|) into @ with R = rl, we obtain a GW Dirac operator 



D 



a r CtC(J + a 2 r 2 C^C)- 1 -C^(l + a 2 r 2 CC^y x 
C(l + a 2 r 2 C^C)- 1 a r CC\\ + a 2 r 2 CC^y l 



(2f 
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where 



C = {<r%) - Wfat^W . (29) 

The locality of D depends on the gauge configuration as well as the value of 
r. If r is zero, then D is equal to D c which is nonlocal, since D c is chirally 
symmetric, free of species doublings and has correct continuum behavior. For 
sufficiently smooth gauge configurations, D is local ( and not highly peaked 
in diagonal elements ) for r within a proper range. For example, consider the 
gauge configuration in Fig. 1, D is local for r e (0.2, 0.8). A rigorous proof of 
the locality of D for a given r and for gauge configurations satisfying a certain 
bound is beyond the scope of the present paper. 
The axial anomaly of D can be written as 

A D (x) = tr[(I + aVCC*)- 1 - (I + aVC+C)- 1 ]^, x) (30) 



where tr denotes the trace over the color and spinor space. Since D c (p6|) in 
the free fermion limit is free of species doubling and has the correct continuum 
behavior, the perturbation calculation in ref. [17] showed that D (28) has the 



correct axial anomaly in the trivial sector. This has been verified explicitly 
on finite lattices. An example is shown in Fig. 1, in which the axial anomaly 
Ad(x) is plotted for each site on a 12 x 12 lattice with lattice spacing a — 1, 
comparing with the topological charge density p(x) = ^Fi 2 (x) of the trivial 
gauge background. The position of a site with coordinates (xi,X2) is repre- 
sented by an integer x = 12 (x 2 — 1) + x±, as the x-coordinate in Fig. 1. The 
axial anomaly Ad(%) is denoted by diamonds, while the topological charge 
density p(x) of the gauge background by circles. The line segments between 
circles are inserted only for the visual purpose. Evidently Ad(x) agress with 
the topological charge density p(x) at each site. 

It is instructive to compare the axial anomaly of the GW Dirac operator 
D (|28f) in a trivial gauge background, as shown in Fig. 1, to that of the 
Wilson-Dirac operator Dw (PS) in the same gauge background, as shown in 



Fig. 1 in ref. The former agrees with the topological charge density 

p(x) while the latter does not. This demonstrates that the transformation 
T(r + R c ) indeed plays an important role in converting Dw into a GW Dirac 
operator D = T{r + R c )[Dy/} which is free of 0(a) lattice artifacts. Thus D 
can reproduce the correct axial anomaly even on a finite lattice. 

However, this GW Dirac operator (p8|) does not possess exact zero modes in 
any topologically nontrivial gauge background ( i.e., D is topologically trivial 
). This is expected since the Wilson-Dirac operator Dw is topologically trivial, 
hence the number of zero modes remains zero ( n±(D) = n±(Dw) = ) under 
the transformation Dw — > D — T(r + R c )[Dw}- Consequently, the sum of 
the axial anomaly over all sites is zero, not equal to the topological charge of 
the gauge background. Therefore the axial anomaly of D ( or Dw ) cannot 
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agree with the topological charge density of the nontrivial gauge background 
on any finite lattice. It follows that the disagreement must persist even in the 
continuum limit a — > W2\. 



Even though the new GW Dirac operator (|28|) is topologically trivial, it 
has well-defined chiral properties and is not plaqued by the notorious problems 
which occur to the standard Wilson-Dirac operator ( p0|) . Therefore, from 
this viewpoint, it seems that (p28[) may be an alternative to the Wilson-Dirac 
operator in lattice QCD, especially for studies involving observables which are 
dominated by the trivial sector. Further implementations of ( f28f) for lattice 
QCD are beyond the scope of this paper. 

For massive Dirac fermion operator, the mass term m should enter D ac- 



cording to the following transformation [IB 



T m (R): D c ^ D= (D c + m)(l + aRD c )- 1 . (31) 

Then for the massive Wilson-Dirac fermion operator, the transformed GW 
Dirac operator with R = rl is 



D 



{ml + ar &C){1 + a 2 r 2 0C)- x -{I- ma r)C\l + aVCC*)" 
{l-ma r)C(I + a 2 r 2 C^C)- 1 {ml + ar C&)(1 + a 2 r 2 CC\ 



(32) 



where C is defined in ( p9|) . 

Although the Wilson-Dirac fermion operator D\y is topologically trivial, it 
can be used to construct the Domain- Wall fermion operator on a five dimen- 
sional lattice with fermion mass as a step function in the fifth dimension such 
that chiral fermions can be realized on the four dimensional walls where the 
mass defects locate. Then the current flowing into the fifth dimension can in- 
duce the correct axial anomaly on the four dimensional walls. Therefore exact 
zero modes can be reproduced for the chiral fermions on the four dimensional 
walls for topologically nontrivial gauge backgrounds, and the Atiyah-Singer 
index theorem can be satisfied. However, the chiral symmetry of the fermions 
residing on the four- dimensional walls is an exact symmetry only in the limit 
the number of sites in the fifth dimension N s goes to infinity. Thus, at finite 
N s , the Domain- Wall fermion may suffer from anomalous effects due to chi- 



ral symmetry violations It seems that one may use the new GW Dirac 
operators (p8 ) and (0) to construct a corresponding Domain- Wall fermion op- 
erator. Then it would possess the exact chiral symmetry ([I]) even at finite N s , 
and the absence of additive quark mass renormalization can be guaranteed. 



4 Summary and discussions 

We can understand the topologically invariant transformation T(r + R c ) by 
the following considerations. 
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Given any lattice Dirac operator D, there exists an operator R such that 
the Ginsparg- Wilson relation ([J) can be satisfied. Then, according to 
there exists a chirally symmetric D c such that D can be written in the form 

D = D c (l + aRD c )- 1 = T(R)[D C ] . (33) 

Thus the inverse transformation of ([£J) is 

D c = T(-R) [D] = D(l — aRD)- 1 = T(R C ) [D] (34) 

where 

R c = -R= + 75^75) • (35) 

Once D c is obtained, it can be substituted into (|3|) with R independent of D. 
In particular, for R = rl, it gives 

D' = T(r) [D c ] = T (r) o T(R C ) [D] = T(r + R c ) [D] . (36) 

If D is topologically proper and satisfies (iii), (iv) and (|), then these 
properties are preserved under the transformation (0), i.e., D' is topologically 
proper and satisfies (iii)-(iv), ([l]) and @. Then it follows that the axial 
anomaly of D' would agree with the topological charge density of the gauge 
background, provided that D' is local ( with r in the proper range ). On the 
other hand, if D is topologically trivial, and satisfies (iii)-(iv) and (|2|), then 
the axial anomaly of a local D' ( with r in the proper range ) would agree with 
the topological charge density for a trivial gauge background, but not for the 
nontrivial ones. 

For the Wilson-Dirac lattice fermion operator (0), the transformed GW 



Dirac operator (|28| ) is local, free of 0(a) lattice artifacts, has correct axial 
anomaly in the trivial sector, and is not plagued by the notorious problems ( 
e.g., additive mass renormalization, etc. ) which occur to the Wilson-Dirac 
operator. It seems that ( p8j) or (^) may be an alternative to the Wilson-Dirac 
operator in lattice QCD. 
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Figure 1: The axial anomaly Ad(%) [ Eq. (|30"D ] of the massless GW Dirac 
operator D [ Eq. flU) ] 

in a trivial gauge background on a 12 x 12 lattice. The 
value of r has been set to 0.5, and there is no significant changes to Ad(x) 
for r e (0.2, 0.8). The background gauge field is the same as that in Fig. 1 of 
ref. []E|. The axial anomaly Ad(x) is denoted by diamonds. The topological 
charge density 5^-^12 °f the gauge background is plotted as circles, which are 
joined by line segments for the visual purpose. 



12 



